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Abstract 

In a recent study we have reported a new type of trial wave function symmetric under the ex- 
change of particles and which is able to describe a supersolid phase. In this work, we use the 
diffusion Monte Carlo method and this model wave function to study the properties of solid 4 He in 
two- and quasi two-dimensional geometries. In the purely two-dimensional case, we obtain results 
for the total ground-state energy and freezing and melting densities which are in good agreement 
with previous exact Monte Carlo calculations performed with a slightly different interatomic poten- 
tial model. We calculate the value of the zero-temperature superfluid fraction p s /p of 2D solid 4 He 
and find that it is negligible in all the considered cases, similarly to what is obtained in the perfect 
(free of defects) three-dimensional crystal using the same computational approach. Interestingly, 
by allowing the atoms to move locally in the perpendicular direction to the plane where they are 
confined to zero-point oscillations (quasi two-dimensional crystal) we observe the emergence of a 
finite superfluid density that coexists with the periodicity of the system. 

PACS numbers: 67.80.-s,02.70.Ss,67.40.-w 
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I. INTRODUCTION 



Quantum crystals are characterized by unusually large atomic kinetic energy, Lindemann 
ratio and non-negligible anharmonicity even at low temperatures and high pressures.-^ The 
counterintuitive possibility of simultaneous solid order and superfluidity in solid 4 He, the 
most representative of quantum crystals, has long attracted the interest of both theoreticians 
and experimentalists. After several unfruitful experiments to detect superfluid signals in 
solid 4 He, Kim and Chan reported few years ago the first evidence of non-classical rotational 
inertia (NCRI) both in confined environment 4 - and in bulk.- From then on, several other 
experimental groups (up to 5, so far) have observed NCRI using different samples containing 
small or ultra small 3 He concentrations, in a simple crystal or in a polycrystal, and using 
several annealing schemes.- There is almost overall agreement of all the data concerning the 
onset temperature T = 75—150 mK at which the superfluid fraction becomes zero, the lowest 
value corresponding to ultra pure samples (only 1 ppb 3 He). However, the experimental 
values of the superfluid density reported so far change by more than one order of magnitude 
(p s /p ~ 0.03 — 0.5%) depending on the purity, annealing conditions in which the crystal is 
grown, etc. Such high dispersion suggests that the superfluid signal observed in solid 4 He 
is probably due to the presence of some defects in the crystal, which could be of different 
nature: dislocations, vacancies or grain boundaries. 7 

On the theoretical side, path integral Monte Carlo (PIMC) calculations performed at 
low temperatures (down to 0.1 K) show that perfect commensurate 4 He crystal does posses 
neither finite superfluid fraction^ nor condensate fraction.- Only non-zero p s /p has been 
estimated in the presence of disorder introduced in the form of a glassy phased and of 
defects like dislocation lines,— and vacancies.— Moreover, our recent calculations based on 
the diffusion Monte Carlo (DMC) method show negligible superfluid fraction of perfect bulk 
solid 4 He at strictly zero temperature.— These DMC calculations have been performed using 
a new model of trial wave function which allows simultaneously for spatial solid order and 
Bose-Einstein symmetry and with the benefit of a simple use for importance sampling. It 
has been shown that the energetic and structural properties of solid 4 He can be reproduced 
very accurately with this trial wave function model when it is used for importance sampling 
in the diffusion Monte Carlo method. 

In this work, we extend our study of solid helium to purely two- and quasi two-dimensional 
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geometries relying upon similar computational approaches to the ones used in Refs. [12IJ13 ]. 
The motivation for carrying out the present study is four- fold. First, relevance of quantum 
fluctuations is generally enhanced in systems of reduced dimensionality hence possible sig- 
natures of superfluidity may be detected more easily. Second, from a computational point 
of view low- dimensional systems are reasonably affordable so one can explore wide thermo- 
dynamic ranges on them. Third, in a recent study^ of strictly two-dimensional solid H2 we 
have shown that when the density of particles is reduced down to practically the spinodal 
point a finite superfluid fraction is observed to appear in the film; in this work we carry out 
similar investigations on solid 4 He (that is, in the metastable regime) in order to unravel 
possible connections between the density of particles and superfluid fraction. And fourth, 
theoretical predictions on low-dimensional model systems can provide valuable understand- 
ing on experimental realization of solid helium confined to restricted geometries and also on 
interpretation of supersolid signatures in general.— - — 

There is previous work done on the estimation of the ground-state properties of strictly 
two-dimensional solid 4 He. Many years ago, Whitlock et al.— performed a systematic study 
of the energetic and structural properties of this system based on the Green's function Monte 
Carlo approach (GFMC). Posteriorly, Gordillo et al.— estimated the phase diagram of two- 
dimensional 4 He over a range of temperatures and coverages using PIMC calculations. More 
recently, various authors have reported on the melting transition and dynamical properties 
of helium films using variational and exact ground-state methods.— 1 ^ Interestingly, Vitali 
et al. 21 have investigated the existence of off-diagonal long range order (ODLRO) in 2D 
solid 4 He using the zero-temperature version of PIMC adapted to the shadow wave function 
formalism, the so-called shadow path integral ground-state method (SPIGS). In the present 
article, we provide comparison with respect to the results reported in these previous works 
and present new predictions as well. In particular, we report on direct estimations of the 
superfluid fraction in 2D solid helium and its dependence on the density of particles. We 
also analyse the superfluid behavior of a quasi two-dimensional crystal, i. e. an ensemble of 
He atoms confined within a plane but allowed to explore the out-of-plane direction locally, 
and assess its dependence on the density of particles and degree of confinement. In the quasi 
two-dimensional case, we observe the presence of a superfluid signal which coexists with the 
periodicity of the system. 

The remaninder of the article is organized as follows. In Section [III we summarize the 
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basics of the DMC method and describe the symmetrized trial wave function model used 
throughout this work. Next, we report results for the ground-state properties of two- and 
quasi two-dimensional solid helium. Finally, we present some discussion and the conclusions 
in Section IIVI 

II. METHOD AND TRIAL WAVE FUNCTION 

We study the ground-state of two-dimensional solid 4 He by means of the DMC method 
and Hamiltonian H = -h 2 /2m He Y^i v ? + T,i<j v ( r a) > with N 

being the number of par- 
ticles. The He-He atomic interaction is modelled with the semi-empirical pairwise potential 
due to Aziz et alM (heretofore referred to as Aziz II). The DMC method solves stochas- 
tically the imaginary-time (r) Schrodinger equation providing essentially exact results for 
the ground-state energy and diagonal properties of bosonic systems within controllable sta- 
tistical errors. For r — > <x>, sets of configurations (walkers) Ri = {ri, ...,tn} generated 
with DMC render the probability distribution function (^f ^), where \l/o an d ^ are the 
ground-state wave function and trial wave function for importance sampling, respectively. 
The short-time Green's function approximation that we use, and according to which the 
walkers evolve, is accurate up to order (At) 3 ; technical parameters in the calculations, as 
for instance the mean population of walkers (= 400) and time step At (= 5 • 10~ 4 K _1 ), 
have been adjusted in order to eliminate possible bias in the total energy per particle to less 
than 0.02 K/atom.-^ 

Customarily, structural and energetic properties of solid 4 He are explored with the 
Nosanow-Jastrow trial wave function 

N N 

* N j(ri, , ix) = Y[ /(r fi ) Yl g(m) = ipjipL , (1) 

i<j i,I=l 

where N is the number of particles (in this work we consider commensurate crystals only 
so iV is also equal to the number of lattice sites), f(r) being a two-body correlation factor 
accounting for atomic correlations, and g(r) a one-body localization factor which accounts 
for the periodicity in the system by linking every particle to a particular lattice site of a 
perfect crystal structure. The wave function \I/ N j leads to an excellent description of the 
equation of state and structural properties of quantum solids^ but it can not be used to 
estimate properties which are directly related to the quantum statistics. The reason of this 
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is that ^nj is not symmetric upon the exchange of particles and it misses the quantum 
statistics of the system. 

In a recent work,— we have introduced a new type of wave function, ^snJj which re- 
produces crystalline order and fulfills Bose-Einstein symmetry requirements simultaneously. 
This model wave function is expressed as 

N N / N \ 

*SNj(ri,...,r^) = n/K)n ' (2) 

i<j J=l \i=l / 

where the product in the second term runs over lattice site indexes. Compact and managable 
analytical expressions for the drift velocity and kinetic energy derive from Eq. ([2]), so \1/snj 
is very well-suited for implementation in DMC codes. This model has proved to perform 
excellently in the description of bulk solid 4 He and also p-H 2 in two dimensions.^ The key 
point in \1/snj is that the localization factor (second term in Eq. |2]) is constructed in such a 
way that voids originated by multiple occupancy of a same site are penalized (this feature 
will be illustrated in brief by a simple example). 

A similar model, \1/lnj, has been proposed recently by Zhai and Wu, 26 which is 

N N / N \ 

* LNJ ( ri , , v N ) = n f( rij ) n 5>(r„) (3) 

i<j i=l \J=1 / 

and where the product in the second term runs over particle indexes. This wave function 
also fulfills quantum symmetry requirements and is well-suited for DMC purposes, however 
it does not account for accurate description of quantum solids. We observe that when 
^lnj is used for importance sampling solid order is not preserved but instead glassy-like 
configurations are generated in the simulations.— In fact, substantially better variational 
energies are obtained in two-dimensional solid hydrogen when using trial wave function 
\&snj instead of \1/lnj (see Table I in Reference 13j). Similar variational outcomes are also 
found in two-dimensional solid 4 He. For instance, at density p = 0.525 a- 2 (a = 2.556 A) 
we obtain £'/iV SNJ = 2.64 (4) K (b = 1.1 a and a = 7.5 ct~ 2 , see next section) whereas 
^yjyLNj _ ^ K (6 = 1.3 cr, a = 7.5 cr -2 and c = 4.0, c appearing in the exponent of 

the McMillan factor, see next section). 

The poor variational quality of wave function ^lnj can be understood in terms of the 
localization factor which, contrarily to what is required to keep solid order, does not penalizes 
multiple occupancy of a same site. By multiple occupancy of a same site here we mean large 
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probability of two particles near a same site to get too close one to the other (that is, as 
it would be allowed by the Jastrow factor alone). Differences between trial wave function 
^snj and ^lnj can be illustrated by a simple example of two particles in a one-dimensional 
lattice. For the sake of simplicity, we assume the distance between the atomic equilibrium 
positions to be one, the parameter in the Gaussian factors (g(r)) a = 1/2 (in arbitrary 
units) and switch off the Jastrow factor. The value of the square of wave function \1/snj 
and ^lnj, |^soz| 2 , in the case of pinning one of the particles in one lattice site (at x — 0) 
and then move the second particle towards it, is plotted in Fig. [TJ As one observes in 
there, the value of \I/lnj at point x — 1 and (which correspond to particles placed over 
different sites and particles placed over the same position, respectively) is identical, whereas 
^snj(^ = 1) > ^SNj(a ; = 0) . Moreover, in the event of atomic overlap (x = 0) the drift 
force 1/^ (d^f/dx) corresponding to wave function ^snj is much more repulsive than that 
of wave function ^lnj- In fact, the curve obtained in the ^lnj case resembles that of a 
liquid where the localizing factor can be thought of a constant. Also it must be noted that 
the value of \I/lnj is maximum at half the way between and 1, thus it will promote larger 
diffusion of the atoms throughout the volume. 

A symmetrized trial wave function that has been successfully applied to the study of solid 
4 He is the called shadow wave function (SWF), proposed by Reatto et al. more than twenty 
years ago.— ^ In the SWF formalism, an array of subsidiary particles (shadow particles) 
is defined and made to interact with the real atoms of the system; shadow particles are 
correlated among them and their coordinates are integrated over the whole volume in such 
a way that bosonic symmetry requirements are fulfilled by construction. At the variational 
level the SWF has been shown to provide very accurate description of solid and liquid helium. 
Nevertheless, this kind of trial wave function has never been used for importance sampling 
in a DMC calculation. In spite of this, very recently the SWF has been implemented within 
the path integral ground state (PIGS) formalism so that variational constraints in principle 
have been removed. This formalism has been used to explore solid 4 He in two dimensions^ 1 - 
and we will comment on their results in the next section. 
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FIG. 1: |^lnj| 2 and I^snj) 2 (Jastrow factor equal to unity) functions in the simple case of two 
particles moving in one dimension and with lattices sites separated by one arbitrary unity. 

III. RESULTS 
A. 2D solid 4 He 

DMC simulations of 2D solid 4 He in the triangular lattice configuration have been carried 
out for N = 120 particles in a (x, y) box where periodic boundary conditions are applied. 
Correlation functions in Eq. ([2]) are chosen of McMillan, f(r) = exp [—1/2 {b/r) 5 } , and 
Gaussian, g(r) = exp [—1/2 (ar 2 )], form. Parameter b and a in factors f(r) and g(r) have 
been optimized using variational Monte Carlo (VMC) at a density of p = 0.480 <j~ 2 (a = 
2.556 A), and we obtain b — 1.1 a and a = 7.5 cr~ 2 as best values (we neglect their 
weak dependence on density). Size effects have been corrected by assuming that atoms 
distribute uniformly beyond half the length of the simulation box. Wave function \1/snj 
already provides a good description of two-dimensional solid 4 He at the variational level; 
for instance, at density p = 0.550 a~ 2 we obtain a variational total energy per atom E/N 
of 3.29 (3) K which must be confronted to the SWF result 2 ^ and variational benchmark 
E/N SWF = 3.10 (1) K, and the Nosanow- Jastrow result E/N NJ = 3.18 (2) K. 
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In Figure [2J, we plot DMC results for the total energy per particle E/N as a function of 
density in solid and liquid 4 He. Results for the liquid phase are taken from Reference |2Qf | - 




Previous Green's function Monte Carlo (GFMC) estimations obtained with a slightly dif- 
ferent atomic pairwise interaction than used here (referred to as Aziz 1)^ are included in 
the plot for comparison. 18 In the solid phase, we fit our total energy per atom results to the 
polynomial curve E/N = e + ap + bp 2 + cp 3 with e = -9.28(0.6) K, a = 87.41(3.4) Ka 2 , 
b = —261.87(6.1) K<7 4 and c = 258.91(3.7) Ka 6 the set of parameters which best reproduces 
them (statistical uncertainties of the fit are expressed within the parentheses). Once the 
energy function E(p) is known in the liquid and solid phases, the corresponding melting 
and freezing densities, namely pi and p s , can be estimated by means of the double-tangent 
Maxwell construction. As a result, we obtain pi = 0.492 a~ 2 and p s = 0.456 a -2 which lie 
in between previous GFMC±£ (pf FMC = 0.471 a" 2 , pf FMC = 0.443 (x~ 2 ) and variational 
SWF-22 (pfWF = Q.522 a' 2 , pf^ F = 0.475 a' 2 ) estimations. The small discrepancy with 
respect to the GFMC results can be understood in terms of the small differences in the 
atomic pairwise potential used. For instance, it is well-known that the Aziz II potential 
provides atomic total energy values around 0.1 K smaller than the Aziz I potential does.— 

A well-known drawback of the NJ model (Eq. [1]) is the impossibility of answering the 
fundamental question whether off-diagonal long range order (ODLRO) and/or superfluid 
behavior may manifest or not in quantum solids. The SNJ model (Eq. [2]) correctly fulfills 
the Bose- Einstein statistics and provides that information to some extent. Quantitatively, 
ODLRO is measured by the condensate fraction no, which is estimated through the asymp- 
totic behavior of the one-body density matrix p(r)/p, namely uq = lim^oo p(r)/p. The 
one-body matrix is an operator which is non-diagonal in coordinate space and does not 
commute with the Hamiltonian H ([H, p] ^ 0) so that DMC output for uq is a mixed esti- 
mator, though bias stemming from the trial wave function can be reduced significantly by 
means of extrapolated estimator techniques.—" 1 ^ 

Vitali et al— have recently adapted the PIGS formalism to the symmetrized SWF to 
investigate strictly two-dimensional solid 4 He with it. Interestingly, the authors of this work 
conclude with the non-existence of ODLRO in perfect 2D solid 4 He. We will comment on 
this finding in the next paragraph in the light of our p s j p results. 

Differently to the estimation of no, the superfluid density of a bosonic system can be cal- 
culated exactly using DMC (whereas this has not been possible yet within the PIGS method) 
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FIG. 2: Energy per particle E/N of solid (S) and liquid (L) 4 He in two dimensions and zero 
temperature. Previous GFMC calculations 18 obtained with the Aziz I pair potential are shown for 
comparison (triangles). Results for the liquid phase are from Ref. [291 ] . 

by extending the winding-number technique, originally developed for PIMC calculations, to 
zero temperature.— Specifically, the expression of the superffuid fraction reads 

El = ^ ma (£M) i (4) 

P T->00 V T / 

where a = N/4D (two-dimensional case) with D = h 2 /2m, D s {r) = ((Rcm(t)— R-cm(O)) 2 ) 
and Rcm is the center of mass of the particles in the plane. In Figure we plot the function 
D s (t) calculated in the solid film at three different densities located near, above and below 
the corresponding freezing density. According to Eq. 01]), the superffuid fraction p s j p can be 
estimated directly from the slope of D s (t) at large imaginary time. In all the studied cases 
we find that the superffuid fraction of perfect two-dimensional solid 4 He is vanishingly small, 
or to be more exact, it lies below our numerical threshold of ~ 10~ 5 . We found analogous 
p s j ' p results in the perfect three-dimensional case. 12 As noted before, Vitali et alr^ have 
recently studied perfect 2D solid 4 He and concluded with the non-existence of ODLRO on 
this system. 

Very recently, we have studied strictly 2D solid H 2 at zero-temperature with analogous 



approaches to the one used here.— In two-dimensional molecular hydrogen, we found that 
in the regime of very low densities (negative pressure regime, pn 2 < 0.390 cr~ 2 ) a finite 
superfluid fraction appears in the crystal. Motivated by this result, and to understand the 
relation between normal and superfluid densities, we have carried out analogous calculations 
in 2D solid helium at stable and metastable conditions (that is, at densities above and below 
pi = 0.492 cr~ 2 , respectively). It is worth noticing that the DMC method has already been 
used to study ground-state properties of metastable liquid 4 He (overpressurized) .— In the 
present case, we find a null value of p a j p for any density down to p = 0.390 er~ 2 (see Figure|3]). 
This result seems to be at odds with our previous findings in H 2 since solid helium possesses 
larger degree of quantumness compared to solid hydrogen. A possible explanation for this 
is that the density p = 0.390 o~ 2 is still far from the spinodal point of two-dimensional solid 
helium. In fact, 4 He is more compressible than H2 (that is, | dP/dV \n e <\ dP/dV |h 2 ) so it 
is likely that the critical density at which mechanical instabilities appear in the first system 
is below that of the second. Also it must be noted that very dilute solid helium films are 
far from being realizable since the liquid phase is always energetically more favourable at 
densities below p = 0.480 a -2 (contrarily to what occurs in two-dimensional H2). In order 
to complete our study of low-dimensional solid helium we have explored a system which can 
be considered as somewhat more realistic, namely a quasi two-dimensional film. 

B. Quasi 2D solid 4 He 

Very recent torsional oscillator-like experiments performed in the second layer of solid 
4 He adsorbed on graphite seem to point towards the possible existence of a new kind of 
supersolid phase. 33 Physical quantities such as the density of particles, temperature and 
degree of corrugation with the substrate, appear to have an important effect on the value 
of this supersolid-like signal. Aimed at investigating on the origins of this manifested low- 
dimensional supersolidity, which is totally absent in strictly two-dimensional solid 4 He, we 
have studied an interesting kind of simple model : a quasi 2D film. In this work, a quasi 
two-dimensional solid refers to a system of interacting 4 He particles with atomic displace- 
ments mostly confined to a plane but which can spread over the z-axis due to zero-point 
oscillations (see Figure H]). This model grasps some essential features of helium layers ad- 
sorbed on carbon-based surfaces like graphene and graphite and could be relevant to describe 
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FIG. 3: Diffusion of the center of mass of two-dimensional solid 4 He in imaginary time calculated 
at a series of densities near, below and above the melting point. A small upwards shift has been 
applied to the curves calculated at higher densities in order to appreciate details of their slope. 

surface effects in bulk crystals. There are several advantages in exploring this model system 
instead of performing more realistic simulations of helium films adsorbed on carbon-based 
surfaces.— ~— First, is the reduction of computational cost which derives from ignoring ex- 
plicit interactions with the substrate and that allows us to explore the superffuid properties 
of the model upon a wide range of conditions. Second, confinement in the z-direction can 
be tuned at wish in order to explore the relation between the superffuid fraction p s /p and 
magnitude of the spatial out-plane fluctuations (Az 2 ) = ((z — (z)) 2 ), which are related to 
the strength of the helium film interactions with the substrate. Certainly, realistic simula- 
tions of 4 He films are necessary to fully understand the origins of supersolid manifestations, 
however, here we assume simplified interatomic interactions and atomic structure (only the 
triangular lattice is considered) in exchange for analyzing possible effects deriving from the 
density of particles and strength of the film-substrate interactions. The density range in 
which we concentrate corresponds to that of non-complete first layer of solid helium ad- 
sorbed on graphite (that is 0.416 < p < 0.745 atom/cx 3 }^, so effects relating to promotion 
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FIG. 4: Top and front views of strictly two-dimensional (left) and quasi two-dimensional (right) 
solid 4 He at density p = 0.515 a~ 2 and zero-temperature. 

of atoms to second or higher layers are not considered. 

Local fluctuations of the atomic positions in the z-direction are achieved by imposing an 
external out-of-plane harmonic potential trap V c (z) oc (z — Zq) 2 , where Zq is the equilibrium 
position of the film in that direction. The Hamiltonian describing the quasi 2D system then 
is H = T + VazizII + V c . The symmetrized trial wave function that we use to describe the 
quasi two-dimensional system is 

N N / N \ N / 1 \ 

*&F(r 1 ,...,rj r ) = n/(r«) II E^DO II ex P (5) 

i<j J=l \i=l / i=l V / 

where r^J is the projection of the vector ri — Rj over the xy-plane, lattice vectors being 
{Rj = (a, b, 0)}, and the value of parameter x is explicitly related to the value of the atomic 
spatial ^-fluctuation by (Az 2 ) = 1/2% . Trial wave function ^snj D ^ s equivalent to \1/snj m 
Eq. 02]) but with additional Gaussian localizing factors on the z-direction; these localizing 
functions correspond to the exact Schrodinger equation solution of a particle moving under 
the action of the harmonic potential field V c . Since the computational technique used in 
this section is the DMC method and the magnitude of the atomic z-spatial fluctuations 
analyzed is fairly small, we have not attempted to construct explicit z-pairwise correlations 
on the trial wave function. It must be noted that these correlations are implicitly taken 
into account by the Jastrow factor contained in ^snj D - Variational parameters contained in 
expression (JSJ) are set to the same value than used in the study of strictly two-dimensional 
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FIG. 5: Diffusion of the center of mass of quasi two-dimensional solid 4 He calculated at the density 
p = 0.470 a~ 2 and different values of the atomic z-spatial fluctuations (Az 2 ). 
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TABLE I: Calculated superfluid fraction (expressed in %) of quasi-2D solid 4 He as a function of 
the density of particles and atomic spatial z- fluctuation (Az 2 ) . 

solid 4 He, made the exception of x which is varied according to the value of the spring 
constant corresponding to the harmonic trap. 

Next, we comment on the superfluid fraction results obtained at several densities and 
^-confinement conditions. In Table I, we report the dependence of the estimated superfluid 
density fraction on the density of particles and the ^-coordinate fluctuation (Az 2 ); in Figure [5] 
we also show the evolution of function D s (t) on imaginary time at density p = 0.470 a~ 2 
and several atomic ^-confinements. It is found that in the (Az 2 ) 1 / 2 < 0.08 a cases (not 
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FIG. 6: Radial pair-distribution function calculated on the xy-plane of the quasi 2D film at different 
densities and fixed amplitude of ^-fluctuations (Az 2 ) 1 / 2 = 0.35 a . At density p = 0.600 cr~ 2 , the 
peaks and valleys of the g X y{f) function turn out to be appreciably less pronounced than in the 
other cases considered. 

shown) the superfluid density is always vanishing, which turns out to be consistent to what 
it is found in the strictly two-dimensional case. On the contrary, when confinement on the 
^-direction is shallow, that is case {Az 2 ) 1 / 2 = 0.35 a, the value of p s j p is always large and 
increases as the density of particles is raised (see Table I). Results in the last column of 
Table I appear to show a competition between in-plane and out-plane interactions; as the 
system is compressed, in-plane repulsive interactions become progressively more important 
so that atoms prefer to spread over the z-direction wherein potential confinement is mild and 
they can move more freely. This has the overall effect of enhancing the superfluid response 
of the system. In the {Az 2 } 1 ^ 2 = 0.18 and 0.22 a cases, estimated p s /p trends on density 
are not that monotonic. At (Az 2 ) 1 ^ 2 = 0.18 a, we see that p s /p first increases when the 
density of particles is raised whereas next it diminishes down to zero-value under further 
compression of the system. This behavior can be understood in terms of the imposed z- 
axis confinement and atomic in-plane interactions as well. When the density of particles is 
first increased, atoms minimize their potential energy by spanning over the ^-axis in order 
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FIG. 7: Calculated radial averaged structure factor S xy (k) of the quasi 2D solid film at p = 
0.515 er~ 2 (A; is in units of cr -1 ). Solid lines represent calculations performed with N = 120 
atoms, dashed lines calculations performed with N = 224 atoms and triangles correspond to 
results obtained for a quasi 2D liquid system (N = 120 atoms). In the lower and upper panels, we 
show results obtained for z-confinement (Az 2 ) 1 / 2 = 0.35 a and (Az 2 ) 1 / 2 = 0.22 a, respectively. 

to keep a distance from their neighbors and move as freely as possible. The total space 
available for the atoms then becomes larger so the effective density of the system becomes 
smaller. The value of p s j p consequently increases. However, when density is further raised, 
out-plane atomic excursions are not favourable any more because large displacements along 
the z-direction require too much energy. The value of p s /p then decreases because atomic 
motion is tightly bound to the xy-plane and the effective density of the system becomes 
large. The competition between in- and out-plane interactions as modulated by the density 
of particles is responsible for the enhancement /depletion of the observed superfluid response 
of the system. At z-axis confinement (Az 2 ) 1 / 2 = 0.22 cr, the trend of p s j p on density exhibits 
an intermediate behavior between that found in the 0.18 and 0.35 o cases. 

Regarding the stability of the quasi 2D solid film, we note that in all the studied cases 
crystal-like order has been found as witnessed by (i) marked oscillating shape of the radial 
pair-distribution functions g X y{r) obtained considering the projection of the atomic positions 
on the xj/-plane (see Figure [6]), and (ii) peaked pattern of the corresponding radial averaged 
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structure factors S xy that scale with the number of atoms (see Figure [7]). Nevertheless, in 
the (Az 2 ) 1 / 2 = 0.35 a cases the film is likely to be a kind of glass system since the value 
of the corresponding maximum S xy (k) peaks are small (see Figure [TJ lower panel) and the 
p s /p values obtained are quite large (last column in Table I) in comparison to the results 
obtained upon tighter z- confinement. In fact, large superfluid fractions (~ 10 — 60 %) have 
been estimated in metastable 4 He glass systems using the PIMC method.— Moreover, we 
have calculated the total energy of a quasi 2D liquid system at same density and (Az 2 ) 1 / 2 
conditions than in the quasi 2D solid system, using a Jastrow factor and Gaussian z-localizing 
factors as importance sampling. It is found that in all the (Az 2 ) 1 / 2 = 0.35 a cases the quasi- 
2D solid film is metastable (that is, Ej~ 2D < E^ 2D ) whereas is stable in the rest of (Az 2 ) 1 / 2 
and p cases (that is, E^ 2D < Ej~ 2T> ). For instance, at (Az 2 ) 1 / 2 = 0.35 a and p = 0.515 ct~ 2 
we obtain E^ 2D = 2.71(2) K/atom and K[nj D = 3.66(2) K/atom, while at (Az 2 ) 1 / 2 = 0.18 a 
and same density we obtain E^~ 2T> = 16.29(3) K/atom and E^ 2D = 15.81(2) K/atom. This 
last outcome, in analogy with the 3D case, appears to corroborate the hypothesis that the 
results obtained at z-confinement (Az 2 ) 1//2 = 0.35 a correspond to quasi 2D glass systems. 

It is worth noticing that although particles in the quasi 2D film are allowed to move in 
the z-direction, there is not enough energy to excite the levels of the transverse confinement 
and the system is kinematically two-dimensional. The radial motion is frozen to zero-point 
oscillations and the magnitude of z- fluctuations therefore is related to the length (Az 2 ) 1 / 2 . 
We illustrate this in Figure M where we plot the atomic z-density profile calculated in the 
quasi 2D solid film at p = 0.470 a -2 and (Az 2 ) 1 / 2 = 0.18 a using the pure estimators 
technique,— and compare it with the corresponding normalized Gaussian z-localizing factor 
entering ^snj D ( case X — 16 cr -2 ). 

IV. CONCLUSIONS 

We have studied two- and quasi two-dimensional solid 4 He at zero temperature by means 
of the diffusion Monte Carlo method and using the recently proposed symmetrized trial wave 
function \I/snj as importance sampling. We have estimated the superfluid density fraction 
of two-dimensional solid 4 He at T = and found that is negligible down to a density of 
p = 0.390 cr~ 2 . Importantly, by allowing the atoms to move along the z-axis we observe the 
appearance of a superfluid response that coexists with the crystalline order of the system. 



16 




FIG. 8: Calculated atomic z-density profile (solid line) in the quasi 2D solid system at p = 0.470 a 2 
and (Az 2 )V2 = 0.18 a. For comparison, we plot the corresponding normalized Gaussian z-localizing 
factor (x = 16 <r -2 ) entering trial wave function ^g^j - Distance is in units of a. 

The magnitude of this response is shown to depend on the degree of 2-axis confinement 
and the density of particles. This finding is valuable for the realization and interpretation 
of more realistic simulations of helium layers adsorbed on carbon-based surfaces, where the 
interactions with the substrate must be taken into account accurately in order to make 
rigorous judgements about the existence of superfluidity and/or ODLRO. In view of the 
present results for the quasi two-dimensional solid, it can be suggested that a well-suited 
system where to observe a finite superfiuid signal is the first layer of 4 He on top of graphene 
or graphite^ which stabilizes in a triangular lattice and possesses relatively small density. 
Work to verify this hypothesis is in progress. 
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